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We consider an exact population transition, defined as the probability of finding a state at a final time being 
exactly equal to the probability of another state at the initial time. We prove that, given a Hamiltonian, there 
always exists a complete set of orthogonal states that can be employed as time-zero states for which this exact 
population transition occurs. The result is general: it holds for arbitrary systems, arbitrary pairs of initial and 
final states, and for any time interval. The proposition is illustrated with several analytic models. In particular 
we demonstrate that in some cases, by tuning the control parameters a complete transition might occur, where a 
target state, vacant at t = 0, is fully populated at time r. 

PACS numbers: 03.65.-w, 32.80.Qk 



Introduction. — The central goal of quantum control is the 
transfer of population from an initial state to a final target state 
HI |2] • Within the framework of coherent quantum control, fo- 
cus has been primarily on designing specific laser-based sce- 
narios that achieve this goal (for some bound state examples 
see, e.g., 0, S]), whereas within the framework of opti- 
mal control, focus has been on identifying control fields that 
achieve this goal, both computationally and experimentally. 

Despite the enormous interest in this area there are very few 
analytic control results about realistic systems. These include 
theorems such as that of Huang-Tarn-Clark a theorem by 
Ramakrisnan et al. on the dimensionality of the Lie Algebra 
induced by the interaction between the system and the control 
field and a theorem by Shapiro and Brumer [8], where 
control was shown to depend on the dimensionality of the con- 
trolled subspaces. As a consequence, any proven fundamental 
result adds considerably to the knowledge base (e.g., 19D). In 
this paper we expose a universal feature of quantum dynamics 
that has significant implications for control. The focus here 
is in the proving this dynamical result; future studies will be 
directed to control applications. 

Specifically, consider an initial state |^(0)) that evolves un- 
der a Hamiltonian H to yield the state |$(t)) at time t = r. 
Of interest is the probability P/(0) = |(/|*(0))| 2 of the sys- 
tem being initially in state \I) undergoing a transition with 
probability Pf(t) = | (^^(r)) | 2 to an orthogonal compo- 
nent \F) at time r. We focus on the possibility of an "exact 
quantum transition" between these states defined as 



P F (r) = Pr(0), 



(1) 



i.e., where the probability of observing state \F) at final time 
t equals the probability of observing the state \I) initially. 

In this paper we prove that there always exists, for arbitrary 
evolution operator and for an arbitrary time r, a complete set 
{'I'fc(O)} of orthogonal states that undergo the exact state tran- 
sition © from | J) to |.F). For a given Hamiltonian, the mag- 
nitude of the associated Pp (r) is determined by the choice of 



r, \I) and \F). As examples, we obtain the set {^(O)} for 
some analytical models, and furthermore provide instances of 
significant transfer, defined by Pp(0) *C -Pr(O). 

While this universality might seem surprising, we show be- 
low that it simply stems from unitarity of quantum evolution. 
Based on unitary evolution, the universal existence of exact 
quantum state transmission between different subspaces was 
demonstrated in [ 10], and cyclic quantum evolution in the the- 
ory of geometric phase was established in II 1 211 . Unitarity is 
also at the heart of the no-cloning theorem, which is funda- 
mental to quantum information science. In the present case 
an inclusive theorem in quantum dynamics based on unitarity 
is derived that is expected to be influential in quantum tech- 
nologies. In particular, note that the dynamical principle is 
here established within the same Hilbert space, unlike ifioll . 
giving an approach that is propitious for a broad range of ap- 
plications, e.g., for quantum computing [11], coherent control 
of atomic and molecular processes []2[], and laser control of 
chemical reaction in molecules HI. 

Universality of the exact population transition. — Consider 
an M dimensional system (M can be infinite), spanned by the 
bases {\a}} and described by density matrix p. The equality 
in Eq. ((TJ becomes 



tr[|F) (F\p(r))=m (I\ p(0)}. 



(2) 



We assume that the system is prepared in a pure state ^ (0) so 
that p(0) = |^(0)) (*(0)|. 

Proposition. There always exists a complete orthogonal 
set {^fc(0)} T , which depends on r, such that an exact popu- 
lation transition described by Eqs. (Q~|) or (|2]i takes place if the 
initially prepared state is a member of this set. 

Proof. Assuming that at time t = the state of the system 
is ^(0), the left side of Eq. (|2]) can be written as 

ti[\F){F\p(r)} = tr[\F)(F\U(r)p(0)UHr)} 

= (V(Q)\UHt)\F)(F\U(t)\V(0)) 
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= tr[|7)(7| p'{r% (3) 

where U(t) is the time evolution operator of the system, and 
we have introduced the exchange operator 

£ = \F)(I\ + \I)(F\+S , (4) 

satisfying £ \I) (I\ £ = \F) (F\, with £ Q = E^y.F I") ("I- 
The exchange operator £ swaps the states \F) and |7) while 
keeping other states intact. It is easy to prove that £ 2 = 1. 
We have also defined the auxiliary density matrix p'(r) — 
|*'(t)) (*'(r)|, with *'(t) = W(t)V(0) and W(t) = 
£U(t). This operator behaves similarly to the time-evolution 
operator. It is significant to note that the operator W(t) is 
unitary, satisfying 

W\t)W{t) = U\t)££U{t) = 1. (5) 

As a unitary operator W(t) can be diagonalized to yield a 
complete set of orthonormal eigenvectors {^(O)},- and ex- 
ponential eigenvalues {exp(i0j.)} r . Any vector \I//.(0) in the 
set thus obeys the eigenequation 

I^(r)* fc (O) = exp(z0 fc )^ fe (O). (6) 

Comparing Eq. (O with Eq. (0), we note that if the state of the 
system at time zero \&(0) is one of the ^fe(0)'s in Eq. {§), then 
the equality tr[|7) (7| p(0)] =tr[|F) (F\ p(r)] in Eq. © holds. 
In other words, an exact population transition occurs between 
states 1 7) and \F) independent of the choice of these states 
other than that they are orthogonal. The result is also valid 
for the exchange operator £ = e ia \F) (7| + e tl3 |7) (F\ + 
£o, where a and (3 are real numbers. In this case the unitary 
condition © translates to W^(r)W(r) = W{t)£^£U{t) = 
1, since £t may not be equal to £. 

The above result is a fundamental attribute of quantum dy- 
namics and should serve as a basic building block in quan- 
tum control theory (see also iflQll and 01211 '). Given an ar- 
bitrary Hamiltonian at an arbitrary time t, and an arbitrary 
pair of states |7) and \F), W(t) can be numerically diagonal- 
ized to obtain its eigenvalue spectrum and eigenstates, giving 
the states |7) and \F) between which the "exact transition" 
-Pj(O) = Pf( t ) takes place. This can be readily done for 
small systems, and we illustrate below several simple eigen- 
problems where the spectrum of W(t) can be analytically 
obtained. However, we emphasize that, unlike specific con- 
trol scenarios, this result is universal, arising only from the 
fact that a unitary operator possesses a complete set of or- 
thogonal eigenvectors. Of particular interest in control sce- 
narios are transitions, which we denote as significant, when 
P/(0) > -Pf(O). Ideally, in control scenarios, we seek ex- 
act transitions that are (what we term) complete, i.e. where 
-Pr(O) = 1 and Pp(0) = 0, so that an initial state, fully pop- 
ulated at time zero, transfers its population to a target state at 
time r. Experience gained from individual sample cases be- 
low sheds light on the theorem and will allow one to assess 
future directions for control applications. 



Example: A two level system. — We discuss three variants of 
the two-level-system (TLS) model. In the first static case the 
theorem holds in a trivial way, though there is no actual pop- 
ulation transition. In the second case a weak time dependent 
perturbation leads to a significant population transfer. The last 
example demonstrates that in a delta-kicked TLS a complete 
transition might take place. The unperturbed TLS model is 
described by the states |0) and |1) of energies Eo and E\ re- 
spectively. Taking into account different types of interactions, 
we explore next the transition between |7) = |0) to \F) = 

I. Static two-level system. — Assuming for simplicity that 
E\ = 0, we obtain the eigenstates of W{t), *±(0) = 
^(|0) ± exp(-iE Q T/2) |1)) with eigenvalues e^±; 4>+ = 
—Eqt/2 and = —Eqt/2 + ir. This case is trivial since 
there is no actual transition during the course of time. How- 
ever, the dynamics is still Hamiltonian and hence (fJI is valid. 

77. Two-level system under a time dependent 
perturbation. — Consider next the two states |0) and |1) 
of energies Eq and E\ respectively, on-resonance with 
a periodic perturbation H'(t) = XVcosujt, where A is 
a parameter characterizing the order of the perturbation 
expansion. Setting again \I) = |0), \F) = |1), we obtain the 
approximate eigenstates of W(t) to first order of A, 

*±(°) ^ 7^ lr ^[-(™os0±l)|O) + |l)]; 

V2(l ± ^ cosO) 

where re ie = iX J Q T dse^ Eo ^ E ^ s cos(ws) (1| V |0); r and 6 
are real numbers. If cos 9 > we obtain the following in- 
equality for the initial state , 5+(0) 

KJ|*+(o)>| « -L(i + ^ COS 0) 

> |<F|* + (O))|«-^(l-£cos0). (8) 

Since |(7|* + (0))| > | (F|* + (0)) |, a significant transfer is 
realized here. For cos 9 < the same inequality holds for 
^_(0). In both cases the initial and final probabilities satisfy 

P J (0) = P F (r)«i(l + r|cose|). (9) 

Figure Q] demonstrates an exact population transfer in the 
present model, computed without approximation for the evo- 
lution of the TLS under a harmonic perturbation. Panel (a) 
demonstrates a significant transition, while for a different set 
of parameters panel (b) shows that at specific times (or for 
a designed time dependent field) a complete transition might 
take place, even for weak perturbations. 

The analytic calculation ©-(O exemplifies a "significant 
transition", i.e. where |(7|^ fc (0))| > |(F|\l> fe (0))|. Having 
demonstrated that population transfer can be achieved, we 
further address the question, particularly relevant to control, 
of what is the maximum achievable significance, defined as 
Pf(0) — 7 J /(0). For the TLS case the significance equals 
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FIG. 1: The dynamics of a two-level system under a cosine perturba- 
tion. Ei - E = u) = 1, A = 0.2, {Q\V\l) = (1|V|0) = 1. (a) An 
exact transition with r = 5; (b) A complete transition with r = tt/X; 
Pi(t) (full line), PF(t) (dashed line). The dotted line marks the po- 
sition of r, where an exact transition takes place, Pi(0) = Pf{t). 



(— Tr[o- z p]), which we want to maximize under the conditions 
that (a) the state is pure, and (b) an exact quantum transition is 
achieved at time r. It can be shown that the condition for ex- 
act transition can be rewritten as r 3 = — s ■ r, where we define 
the vectors s(r) = -|Tr [W \t)o z U "(t)ct] and r = Tr[pa]; 
Ci=i,2,3 are the x, y, z Pauli matrices, respectively. If, for ex- 
ample, the component S3 was equal to 1, the condition for an 
exact transition forces r% to be zero, and there is no initial 
state for which there is significant exact transfer. Assuming 
that S3 ^ 1, on the other hand, the significance becomes max- 
imal for the initial state 



*(Q) 



1 
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[(ai-w 2 )|0) + (l-a 3 )|l>] • (10) 



Note then that generally, states that maximize the significance 
will not be eigenstates of W(t). Thus, although in accord 
with the above theorem one achieves exact transitions, the 
ideal complete transition is obtained rarely. However, the 
proven theorem provides a new framework in within which 
to modify the Hamiltonian to achieve transitions with increas- 
ingly larger significance. 

///. Kicked two-level system. — For the same model, again 
with \I) = |0) and \F) = |1), consider a non-perturbative 
time-dependent Hamiltonian, 



H 



w{a z cose + a x sine),0 < t < d 
ujo z ,t > <5 



(11) 



where a x , a y and a z are the Pauli matrices. Here one can 
generically write the operator W(t) = —i(cos9 + ilt ■ 
!? sin 9), where in this case the parameters are 



cosf 
n z sin ( 
n y sin ( 



sin(tn^) sin(e) cos[oj(t — 5)], 
sin(tn^) sin(e) sin[w(r — 5)], 
— cos(nj(5) sin[w(r — 5)] 
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FIG. 2: The population Pi(0) = Pf(j) for a kicked two-level 
model, showing that with proper tuning of parameters a complete 
population transition can be achieved, t — n, e = n/2, ui = 1, 
w = 1. I J) = |0) (full line); \I) = |1> (dashed line), 



— cos(e) cos[w(t — 5)] sin(wS), 
n x sin9 = — cos(e) sin(wS) sin[o;(T — 5)], 

+ cos(wS) cos[w(r — 5)], (12) 

with + riy + n 2 z = 1. The two eigenstates of W(t) are 
*+(0) = ^ 1 ^(-i)e n |0) + |1) with eigenvalue 

5^ |0)- 



and*_(0) 



l-n, , 
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1) with eigen- 



value — ie , where 7 = arctan(n a; /nj,). Hence, adopting 
^+(0) as the time-zero state, we obtain the results 



P f (t) = Pj(0) = 



1 



(13) 



As an example, if e = it/ 2 and a strong pulse a x kicks at 
t = such that zuS — > 7r/2, in the limit 5 — > one gets 
that n z -> 1, *+(0) -> |0) and J7(r)*+(0) -> |1), which 
is a complete quantum transition from |0) to |1), satisfying 
|(J|*(0))| = 1 and |(F|*(0))| = 0. In Fig. E]we show 
that by carefully tuning the interaction parameters, e.g. the 
delay time S, one can achieve such a complete transition. We 
next demonstrate that a complete transition can take place in 
general in an adiabatically evolving system. 

Complete transitions: Adiabatic evolution. — Consider a 
time-dependent Hamiltonian H(t). If it is varied sufficiently 
slowly, the evolution of the system is adiabatic, and the sys- 
tem occupies an (instantaneous) eigenstate of the Hamiltonian 
H(t), provided the time-zero state 1^(0)) is an eigenstate of 
H(0). If the state of the system at time r, |^/(r)), is orthogo- 
nal to the time-zero state, one can obtain a complete transition 
by setting \I) = |*(0)) and |F) = |*(r)), both eigenstates of 
W(t). As an example, consider the magnetic Zeeman effect 
where a magnetic field splits the atomic (or molecular) degen- 
erate levels, characterized by the magnetic quantum numbers 
M. The Hamiltonian is effectively given by 



H = B(e)J z + T(e)J x , 



(14) 



where the second term T(e)J x is responsible to quantum 
transitions between different values of M (— — J, ...J). The 
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time dependent modulation is controlled by the parameter 
e = § — t, and we manipulate the magnetic field such that 
T(e) [5(e)] is an even [odd] function of e, T(±§) = and 
-B(O) > 0. We now choose the initial state as \I) = \ — J), the 
lowest eigenstate of H(0). If we control the evolution of H(t) 
adiabatically from time to r, the state of the system at time 
t becomes \F) = | J), which is the lowest eigenstate of H(t). 
Since | J) and | — J) are orthogonal, the quantum transition 
(fZ|l is complete. The results of Ref. | j 1 3|] may be an example 
of this scenario when J = 1/2. 

Superposition: Eigenstates in a three-level system. — Fi- 
nally, we address the following conceptual question: can one 
achieve an exact population transition ([TJ with a superpo- 
sition of the eigenstates of W{r) as the time-zero state of 
the system? We explain next the conditions for this transfer 
by considering a three-level Hamiltonian H = O ( 1 0) (1| + 
|1) (2| + h.c). If one chooses the initial and final states to 
be \I) = |0) and \F) = |2), the exchange operator £ — 
|0) (2\ + |2) (0| + |1) (1| commutes with H fjj|. One can eas- 
ily obtain the eigenstates and eigenvalues of W{t) = £U(t), 



*«(0) 
*&(0) 

*c(0) 



V2 



(|0)-|2)); e*- 



-(|0)+ 7b |l) + |2»; e*=e 
l(|0)+ 7c |l} + |2)); e ^=e^ n \ (15) 



7h. c are functions of r and ft but their exact form is not im- 
portant for the discussion below. Examining (fTBI l, it is ob- 
vious that there is no significant transition if the time-zero 
state is any of these three eigenstates, since | (/^^(O)) | = 
| (i ;l |4' / r.(0)) | for k = a, b, c. We show next that under some 
strict conditions, a superposition state can yield a complete 
transition. Since the set {^(O)} is a complete orthogonal 
set, a general time-zero state can be expanded as ^(0) = 



Sfcli Cfc^fe(O), not necessarily an eigenstate of W(r), 



M 



W{t)*(0) = ^C fc e^ fc * fe (0). 



(16) 



fe=i 



An exact transition can still take place at a specific time T, 
obeying the eigenvalue equation 



W(T)V{0) = e^ (T) *(0). 



(17) 



This equality is satisfied if C / t[exp(i0/ c ) — exp(i</>(T)] = 0. 
Thus, for the contributing k coefficients, CV- 7^ 0, we obtain 
a set of conditions <f>(T) = (f>k + 2nKk, where Kk are arbi- 
trary integers. This is a very restrictive condition when there 
are many nonzero coefficients. However, it may be still satis- 
fied for particular systems with special symmetry. We exem- 
plify this within the three-level model presented above [See 
Eq. (TT3T>1 . Assuming the following superposition state at time 
zero, <f(0) = C a ^a(0) + C b ^ b (0), we obtain 



W(t)*(0)=- C a * a (0)-e 



C b ^ b (0) 



(18) 



which is generally not an eigenstate of W(r). However, Eq. 
(TTrb is satisfied at the specific time r = ir/ (V2fi) leading to 



(19) 



(/|*(0)) = C b /2-C a /yf2 
(F|*(0)) = C b /2 + C a /V2, 



manifesting a significant transition and a population transfer 



Pj(0) = P F (n/V2n) 



(20) 



\Ca\ 



\C b 



2V2 



(C*C b + C b C a ) 



The transition can be made complete, depending on the super- 
position preparation coefficients C b and C a . 

Conclusion. — We have proved the universality of exact 
quantum transitions, demonstrating that for a given Hamilto- 
nian and a pair of states \I) and \F) in the associated Hilbert 
space, there always exists a complete set of orthogonal states 
that when employed as the time-zero state of the system, lead 
to an exact population transition between the pair |7) and 
\F). This universal proposition is a fundamental feature of 
quantum dynamics and a promising building block in quan- 
tum control. We have demonstrated the result analytically on 
the time-modulated two-level-system model, showing that in 
some cases a complete population transfer can be obtained. 
We have also shown that an adiabatic evolution can lead to a 
complete transition. Finally, we have analyzed exact popula- 
tion transitions in a superposition of states. Applications to 
specific control studies are the subject of future work. 
Acknowledgments. LAW has been supported by the Iker- 
basque Foundation. DS acknowledges support from the Uni- 
versity of Toronto Start-up Fund, and PB was supported by 
NSERC. 



[1] S. A. Rice and M. Zhao, Optical Control of Molecular Dynam- 
ics (Wiley, New York, 2000). 

[2] M. Shapiro and P. Brumer, Principles of the Quantum Control 
of Molecular Processes (Wiley, New York, 2003). 

[3] E.g., P. Krai, I. Thanopulos, and M. Shapiro, Rev. Mod. Phys. 
79, 53 (2007) and E. A. Shapiro, V. Milner, and M. Shapiro, 
Phys. Rev. A 79, 023422 (2009) and references therein. 

[4] E.g., J. B. Gong and S. A. Rice, J. Chem. Phys. 120, 5117 
(2004); J. B. Gong and S. A. Rice, J. Chem. Phys. 120, 9984 
(2004). 

[5] U. Gaubatz, P. Rudecki, S. Schiemann, and K. Bergmann, J. 

Chem. Phys. 92, 5363 (1990); K. Bergmann, H. Theuer, and B. 

W. Shore, Rev. Mod. Phys. 70, 1003 (1998). 
[6] G. M. Huang, T. J. Tarn, and J. W. Clark, J. Math. Phys. 24, 

2608 (1983). 

[7] V. Ramakrishna et al, Phys. Rev. A 51, 960 (1995). 

[8] M. Shapiro and P. Brumer, J. Chem. Phys. 103, 487 (1995). 

[9] For an extension of this result to open quantum systems see, 

L.-A. Wu, A. Bharioke, and P. Brumer, J. Chem. Phys. 129, 

041105 (2008). 

[10] L.-A. Wu, Y.-x. Liu, and F. Nori, Phys. Rev. A 80, 042315 
(2009). 



[11] M. A. Nielsen and I. L. Chuang, Quantum Computation and 
Quantum Information, (Cambridge: Cambridge University 
Press 2000). 

[12] L.-A. Wu, Phys. Rev. A 50, 5317 (1994). 

[13] S. Ashhab, J. R. Johansson, and F. Nori, Phys. Rev. A 74, 



052330 (2006). 

[14] £ does not commute in general with the Hamiltonian. For 
ample, if £ = |0) <1| + |1) <0| + |2) <2|,then [£,H] 0. 



